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ABSTRACT 

We use the ROSAT North Ecliptic Pole (NEP) survey to construct a small, but 
purely X-ray flux- limited sample of cataclysmic variable stars (CVs). The sample 
includes only 4 systems, 2 of which (RX J1715.6+6856 and RX J1831. 7+6511) are new 
discoveries. We present time-resolved spectroscopy of the new CVs and measure orbital 
periods of 1.64±0.02 h and 4.01 ±0.03 h for RX J1715.6+6856 and RX J1831. 7+6511, 
respectively. We also estimate distances for all the CVs in our sample, based mainly 
on their apparent brightness in the infrared. The space density of the CV population 
represented by our small sample is l.lt~ !7 x ^ 5 pc~ 3 . We can also place upper limits 
on the space density of any sub-population of CVs too faint to be included in the NEP 
survey. In particular, we show that if the overall space density of CVs is as high as 
2 x 10~ 4 pc~ 3 (as has been predicted theoretically), the vast majority of CVs must be 



fainter than Lx 



2 x 10 29 ergs 1 . 
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1 INTRODUCTION 

Cataclysmic variable stars (CVs) are semi-detached inter- 
acting binary stars, consisting of a white dwarf primary ac- 
creting from a companion th at is usually a l ate- type, approx- 
imately main-sequence star. IWarnerl (|l995l ) gives a compre- 
hensive review of the subject. 

Mass transfer, and the evolution of CVs, is driven by 
angular momentum loss from the orbit. At long orbital peri- 
ods (P or b ^ 3h), angular momentum loss through magnetic 
braking is thought to dominate over that resulting from 
gravitational radiation. The canonical disrupted magnetic 
braking model of CV evolution was designed to explain the 
period gap (a pronounced drop in the number of CVs at 

2 h < P or b < 3 h). In this scenario, magnetic braking stops 
when P rb reaches ~ 3h, the secondary loses contact with 
its Roche lobe, and gravitational radiation i s left as the only 
angu l ar momentum loss me c hanism (e.g. Robinson et al.l 
1 19811 : IRappaport et all Il983l : ISpruit fc Ritteil 1 19831 ). The 
secondary regains contact with its Roche lobe when grav- 
itational radiation has decreased P or b to ~ 2 h, and mass 
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transfer resumes. As the secondary star loses mass, its ther- 
mal time-scale eventually exceeds the mass-transfer time- 
scale (even though the mass-loss time-scale increases as M2 
decreases). When this happens, the secondary is not able 
to shrink rapidly enough in response to mass loss, so that 
the orbital evolution slowly moves back through longer pe- 
riods (e.g. iPaczyriskil Il98ll : iPaczvhski fc Sienkiewica Il98ll : 
IRappaport et alll 19821 ). CVs in this final phase of evolution, 
where P or t is increasing, are referred to as 'period bounc- 
ers'. This reversal in the direction of change in P or b causes a 
sharp cut-off in the P or b distribution of hydrogen-rich CVs 
at about 76 min, called the period minimum. 

Several serious discrepancies between the predictions of 
theory and the properties of the observed CV sample remain. 
Theoretical models persistently place the period minimum 
about 10 min short o f the observed cut-off at P or b — 76 min 
(see e.g. iKolbl |2002| for a review). It has also long been 
suspected that the predicted ratio of short- to l ong-period 
CVs is larger than is implied by observations (e.g. iPattersorJ 
1998). [Pretorius et al.l (120071 ) recently confirmed that this is 
a real discrepancy, rather than the result of observational 
bias. 

In order to make progress, it is necessary to place quan- 
titative constraints on other key parameters predicted by 
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CV formation and evolution scenarios. The most fundamen- 
tal of these is the overall space density of CVs, p. However, 
observational estimates of p differ by more than an order 
of magnitude at present, and little effort has been made to 
quantify the statistical and systematic uncertainties affect- 
ing these estimates. In order to understand this situation, it 
is useful to take a closer look at some of the more extreme 
valu es in the literature. 

ICieslinski et~ai1 (|2003l ) suggested a very low space den- 
sity of p ^ 5 x 10~ 7 pc~ 3 for dwarf novae (DNe), the most 
common type of CV. This result was based on DN out- 
bursts detected in the OGLE-II survey, and assumes that 
most DNe in the OGLE-II fields within 1 kpc would have 
been recovered. However, the time coverage of the survey 
(~ 4 y) was much shorter than the longest DN recurrence 
time-scales (e.g. ~ 30 y for WZ Sge), and long inter-outburst 
intervals are known to be associated precisely with the in- 
trinsically faintest CVs that probably dominat e the overall 
popul ation. At the other end of the spectrum, IShara et al.l 
(119931 ) obtained a rather high estimate of p ~ 10 _4 pc -3 from 
a photometric survey for faint (possibly hibernating) CVs, 
but here the CV nature of the o bjects that were cou nted 
remains to be established. Finally. [Schwope et al.l l|2002r ) ob- 
tained ~3x 10 -5 pc -3 from a CV sample constructed from 
the ROSAT Bright Survey. However, this number was pri- 
marily based on two systems with estimated distances of 
~ 30 pc, both of wh ich were subsequently sh own to be 5-10 
times more distant l|Thorstensen et al1l2006l ). 

These examples illustrate that the statistical errors 
associated with observational space density estimates are 
usually dominated by uncertain distances and small num- 
ber statistics, whereas the dominant systematic errors are 
caused by selection effects. We also suspect that the hard- 
to-quantify systematic uncertainties usually outweigh the 
statistical errors, which probably explains why virtually no 
space density estimate in the literature has come with an 
error bar. Given all this, our goal here is to provide a new 
determination of p from a small CV sample with very simple, 
well-defined selection criteria. Since our sample does not suf- 
fer from unquantifiable systematic errors, we are also able to 
provide meaningful errors on p that fully account for small 
number statistics and uncertain distances. 

An attribute that should be shared by all CVs is X- 
ray emission generated in the accretion flow. This is useful, 
because it is much easier to construct a reasonably deep flux- 
limited CV sample in X-rays than it is in the optical. For 
an optical sample, it is usually necessary to cut down on the 
number of sources that can be followed up by introducing 
additional selection criteria based on, e.g., colour, variabil- 
ity, or emission lines; any such selection makes it more diffi- 
cult to understand the completeness of the resulting sample. 
Another advantage of an X-ray flux-limited sample is that 
Fx/F pt decreases with increasing P or b (i.e. with increasing 
optical luminosity). This means that CVs do not span as 
wide a range in X-ray as in optical luminosity. 

The ROSAT North Ecliptic Po le (NEP) survey (e.g. 
iGioia et al.1 120031 ; iHenrv et al.l 120061 ) is constructed u sing 
ROSAT All Sky Survey (RASS; see IVoges et"allll999l and 
IVoges et alJl200Ch data around the north ecliptic pole; the 
RASS reaches its greatest sensitivity in this 81 sq.deg. re- 
gion, centred on a 20 oo = 18 h 00 m , <5 2 ooo = +66°33'. The NEP 
survey thus comprises the deepest wide-angle, contiguous re- 



gion ever observed in X-rays; it consists of a complete sample 
of 442 sources above a flux of about 10~ 14 erg cm _2 s _1 in the 
0.5 to 2.0 keV band (the flux limit varies from 1.2 x 10~ 14 to 
9.5 x 10~ 14 erg cm _2 s _1 over the survey are43). This unique 
combination of depth and breadth, together with moder- 
ate Galactic latitude and extinction, provides an excellent 
opportunity to investigate the space density of CVs. Very 
importantly, out of the 442 X-ray sources included in the 
survey, all but 2 (the number of spuri ous sources expected 
in a sample like this; IGioia et al.ll2003T ) have been optically 
identified and spectroscopically observed. Therefore, we are 
able to construct a complete, purely X-ray flux-limited CV 
sample from the NEP survey. 

We use the CVs detected in the ROSAT NEP survey 
to provide robust observational constraints on the CV space 
density, by carefully considering the uncertainties involved. 
The constraints include an estimate based on the observed 
sample of CVs, as well as upper limits on the size of a hy- 
pothetical population of CVs with X-ray luminosities low 
enough to evade the survey entirely. We also present radial 
velocity studies of the two new CVs discovered in this sur- 
vey. 



2 CVS IN THE ROSAT NEP SURVEY 

The NEP survey area includes 3 known CVs, namely IX 
Dra, SDSS J173008.38+624754.7 (we will abbreviate this 
to SDSS J1730), and EX Dra. All three were detected in 
the RASS, but IX Dra is slightly too faint to be formally 
included in the NEP sample. 

We inspected the optical spectra of all Galactic NEP 
sources, and found only two additional objects with the 
spectral appearance of interacting binaries. A more detailed 
study of stellar sources included in the survey has yielded 
accurate spectral c lassifications of all stellar NEP sources 
<|Micela et all 120071 ). and no additional CV candidates. We 
are therefore confident that we have identified every CV in 
the NEP sample. 

The two interacting binary candidates are RX 
J1715. 6+6856 and RX J1831. 7+6511, which were confirmed 
as CVs by the follow-up observations described below. Fig. [T] 
gives finding charts for the two newl y discovered CVs; their 
coordinates are listed in Table [1] Gioia et al.l (120031 ) re- 
port photographic magnitudes of E = 18.6 and O = 18.3 
for RX J1715.6+6856, and E = 14.8 and O = 15.2 for 
RX J1831. 7+6511 (the errors in these measurements are 
±0.5 mag; O and E are similar to Johnson B and Kron- 
Cousins R, respectively). 



2.1 Observations of the newly discovered CVs 

In order to confirm that RX J1715. 6+6856 and RX 
J1831. 7+6511 are CVs, and to measure their orbital periods, 
we obtained medium-resolution, time-resolved spectroscopy 
with the Calar Alto Faint Object Spectrograph (CAFOS) 
on the Calar Alto 2.2-m telescope. Table [T] gives a log of 
the observations. The G-100 grism was used in combination 



Assuming a 10 keV thermal bremsstrahlung spectrum; this 
choice is justified in Section [3] 
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Figure 1. 4' X 4' finding charts of the optical counterparts of RX 
J1715.6+6856 (left) and RX J1831. 7+6511 (right), made using 
white light acquisition images taken with CAFOS. North is at 
the top and east to the left in both images. 



with a slit width of l'.'2, yielding a spectral resolution of 
~ 4.2 A over the wavelength range 4 300 to 8 000 A. A few 
spectra of each object were also taken using the B-200 grism, 
and a slit width of l'.'b. These spectra have lower resolution 
(~ 9.4 A) and were not used for radial velocity measure- 
ments. Th e spectrophotometric standard star BD+28°4211 
l|Okelll990l) was observed on all three nights to provide flux 
calibration. Except for the first hour of the first night, condi- 
tions were photometric during our observations. Regular arc 
lamp exposures were taken to maintain an accurate wave- 
length calibration. 

The data were reduced using stan dard procedu res in 
IRAF0, including optimal extraction (|Hornel Il986l ). We 
computed radial velocities of both the Ha and H/3 lines, 
using the Four i er cr oss correlation method described by 
iTonrv fc Davisl (| 1979), as implemented in the fxcor routine 
in IRAF. In the case of the Ha line, we correlated the wave- 
length range 6 450 to 6 650 A, while for the measurements 
of H/3 radial velocities, the correlation was restricted to the 
range from 4 700 to 5 020 A. For both RX J1715. 6+6856 and 
RX J1831. 7+6511, the Ha line yielded higher S/N radial ve- 
locity curves; we therefore use the radial velocities measured 
for this line in determining the orbital periods. 



2.1.1 Observational results for RX J 1715.6+6856 

Fig. [5] shows the average spectrum of RX J1715. 6+6856, 
which displays broad, double peaked Balmer and He i emis- 
sion lines superimposed on a blue continuum; HenA4686 
and FenA5169 are also detected in emission. The system 
has the spectral appearance of a quiescent DN. No large am- 
plitude brightness variations were evident in our data, and 
we estimate V = 18.3 ± 0.2 from our spectroscopy (where 
the uncertainty is dominated by slit losses). 

At V = 18.3, RX J1715.6+6856 is close to the magni- 
tude limit of our telescope/instrument combination, and the 
time resolution and base line of our data yield only a fairly 
uncertain period. The Fourier transform in Fig.[3]shows that 
the radial velocity curve is aliased; the largest amplitude sig- 
nal is at 0.1695 mHz (corresponding to a period of 1.639 h), 
and is marked by a vertical bar in the lower panel of Fig. [3] 
The two nearest peaks are at 0.1575 and 0.1816 mHz. We 
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Figure 2. The average of all the spectra of RX J1715. 6+6856 
taken with the G-100 grating. Individual spectra were shifted to 
the rest frame using the radial velocities calculated from the Ha 
line before being averaged. The system is identified as a CV by 
its broad, double peaked Balmer, He I, and Hell emission lines. 
Also present is the Fe II A5169 line. 




0.2 0.3 
Frequency/rnHz 

Figure 3. Lower panel: The Fourier transform of the Ha radial 
velocity curve of RX J1715. 6+6856 obtained on 28 and 29 July 
2006. A vertical bar at 0.1695 mHz (1.639 h) indicates the period 
that provides the best fit to our data. Upper panel: The radial 
velocity measurements, folded in the best-fit period. One cycle is 
repeated in this plot; the superimposed sine function is a non- 
linear least squares fit to the data. 



fitted sinusoids to the data by non- linear least squares, using 
the periods and amplitudes of these three signals as start- 
ing values. We found that the 1.639 h period provided a 
significantly better fit than the other aliases (the standard 
deviation of the residuals for a fit near this period is ~ 15 
% lower than for either of the other choices). The 1.639 h 
period also yields the largest best-fit amplitude (66.1 km/s). 
Nevertheless, there is some uncertainty in the cycle count. 
Assuming that we have identified the correct alias, the or- 
bital period is of RX J1715. 6+6856 is 1.64 + 0.02 h. The top 
panel of Fig.|3]shows the radial velocities of Ha, phase-folded 
using this period, with the sine fit over-plotted. 
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Table 1. Log of the observations, and coordinates for the optical counterparts of the two X-ray sources. 



Object 


02000 


^2000 


Date at the start 
of the night 


First integration 
HJD 2453900.0+ 


Grism 


integration 
time/s 


Number of 
spectra 


RX J1715.6+6856 


17:15:41.7 


+68:56:43 


2006 Jul 27 


44.36780950 


B-200 


1000 


7 








2006 Jul 28 


45.37354850 


G-100 


1000 


11 








2006 Jul 29 


46.35562975 


G-100 


800 


9 


RX J1831. 7+6511 


18:31:44.4 


+65:11:32 


2006 Jul 27 


44.51118822 


G-100 


600 


16 










44.66093472 


B-200 


600 


2 








2006 Jul 28 


45.59054418 


G-100 


600 


10 








2006 Jul 29 


46.47792121 


G-100 


600 


22 



2.1.2 Observational results for RX 31831.1+6511 

RX J1831. 7+6511 was descending into a photometric low 
state over the three nights of our observations, fading by 
roughly a magnitude in V. This was most likely the decline 
from a dwarf nova outburst. The spectra show fairly narrow 
Balmer and He i emission lines (the FWHM of the Ha line 
is ~ 800km/s). The H/3 and H7 line profiles change from 
broad absorption wings with emission cores in the bright 
state, to pure emission in the faint state. This transition 
is illustrated in Fig. 2] which displays nightly averages of 
the G-100 spectra after shifting to the rest frame using the 
radial velocities calculated from the Ha line. We also detect 
FenA5169 in RX J1831. 7+6511, increasing in strength as 
the system fades. 

The lower panel of Fig. [5] shows a Fourier transform of 
the radial velocity curve of RX J1831. 7+6511 (including all 
3 nights of data). A vertical bar at 0.0693 mHz (4.01 h) 
indicates the highest amplitude signal, but there is also 
strong power present at the one-day aliases (0.0578 mHz 
and 0.0808 mHz). Non-linear least-squares fits of sine func- 
tions to the data produce the largest amplitude for the 
4.01 h period; the standard deviation of the residuals is also 
roughly 45% larger for fits using the other two candidate 
periods. We are therefore reasonably confident that the sig- 
nal at 4.01 h represents the orbital modulation, and find 
P orb = 4.01 ± 0.03 h. The top panel of Fig. shows the Ha 
radial velocity curve, phase-folded using the orbital period, 
together with the best-fit sine function. The ampliude of the 
fit is 79.7 km/s. 



2.2 Distance estimates 

Several independ ent methods of es timat ing distances to CVs 
are discussed by iPattersonl (|l984h and IWarnerl l|l987l) . We 
will use methods based on DN outburst maximum, the 
strength of emission lines, and the brightness of the sec- 
ondary star. 

The relation between orbital period an d the absolute V 
magnitudes of D Ne at maximum found bv lWarnerl (|l987h . 
as recalibrated bv lHarrison et all {2004!), is 



Ah 



= 5.92 - 0.383P orb /h. 



(1) 



In this equation, Mymax is corrected for binary inclination 
(i) using 

AM v (i) = -2.50 log [(1 + 1.5 cos i) cos i] (2) 
{Paczvnski fc Schwarzenberg-Czernvl Il980l ) . To apply this 
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Figure 4. The average nightly spectrum of RX J1831. 7+6511. 
The system faded from V = 15.0 to V = 15.9 over the three 
nights, with the strength of the emission lines increasing, as is 
commonly seen durin g the decline from dwarf nova outbursts (e.g. 
iHessman et al.lll984h . 



relation one thus needs knowledge of i, as well as a mea- 
surement of V at maximum. The scatter of data around the 
Mvmax-Porb relation is roughly ±0.5 mag. 

The equivalent widths (EW; throughout this paper we 
will take the EWs of emission lines as positive) of disc emis- 
sion lines are theoretically ex pected to increase with decreas- 
ing mass transfer rate (M): [Patterson! (1 19841 ') gives an em- 
pirical relation 

0.55(M v -4) 



EW(H/3) = 0.3M^ 



+ e" 



(3) 



that predicts My.disc (the absolute V magnitude of the 
disc) to within about 1.5 mag for EW(Hf3) > 15 A. For 
EW(H/3)< 15 A, the data imply only My.disc < +6. 
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Figure 5. The low frequency end of the Fourier transform of the 
Ha radial velocity curve of RX J1831. 7+6511 (bottom panel). A 
vertical bar at 0.0693 mHz (4.01 h) marks the signal that most 
likely corresponds to the orbital period. The top panel shows the 
phase-folded radial velocity curve, with the best-fit sine function 
over-plotted (all data are plotted twice). 



iBailevi (| 198ll ) describes a method of estimating dis- 
tances to CVs that relies on finding the -/f-band surface 
brightness of the secondary from V — K colour of the system. 
It has since been shown that this method is not as reliable 
as was first thought, and that the V and K magnitudes of 
a CV alone canno t provide a secu re distance estimate (e.g. 
lBeuermann|[2OO0l ). iKniggel ([2006) provides an alternative, 
based on a semi-empirical donor sequence for CVs that gives 
absol ute magnitud es for the secondaries as a function of 
P orb . IKniggel <|2006h finds offsets of A J = 1.56, AH = 1.34, 
and AK = 1.21 between his sequence and the absolute JHK 
magnitudes of a sample of CVs with parallax distances. The 
scatter of the sample magnitudes around the offset sequence 
are (taj = 1.25, uah = 1.12, and ctak = 1.07. We will use 
these offsets and errors to obtain distance estimates; the 
limitations of the method are discussed in more detail by 
IKniggel iM). 

The infrared observations u sed here are from the Two 
Micron All Sky Survey (2MASS; Iskrutskie et aLll2006r i. We 
interpolate linearly on the high resolution version of the 
donor sequence to find predicted absolute magnitudes for the 
secondaries at the o rbital p e riods of the observed CVs. The 
CIT magnitudes of IKniggel ((2006! ) were transformed to th e 
2MASS system using the prescription of ICarpenterj l|200lf ). 
updated to reflect the final 2MASS data releasqj. Errors 
from the transformation are insignificant for our purposes, 
and will not be quoted (although we did include them). 

To estimate extinction effects, we compared a model 



of Amores fc Lepinei ( 2005^3 to the model of lDrimmel et all 
|2003l ) at several positions in the NEP region, and found 
an average difference of ~45 % be tween the Ay values pre - 
dicted by the two. We will use the I Amores fc Lepinei (|2005l ) 
extinction values, and assume errors of 50 %. Extinction is 
obviously a fuction of distance; we therefore do a few itera- 
tions to ensure that the value we finally use in the distance 
calculation is that given by the model at the estimated dis- 
tance to the object; for lower limits on distances we use 
the total Galactic extinction. To convert from visual extinc- 
tion to that in the 2MASS bands, we use Aj = 0.282Ay- 
A H = 0.175A V -, and A Ks = 0.112Av (jCambresv et all 
I2OO2I) . 

We assume that errors in absolute magnitude, appar- 
ent magnitude, and extinction are Gaussian, and use this 
to find a probability distribution for the distance to each 
source. The distances we quote are in all cases the median, 
together with the symmetric 1-a confidence interval around 
the median (this is the interval defined so that the integral 
of the distribution between its lower bound and the median, 
as well as between the median and its upper bound, is 34% 
of the integral over all distances). 

Two of the three techniques we use to estimate distances 
require the orbital period of the CV. For RX J1715.6+6856, 
RX J1831. 7+6511, and SDSS J1730 there is some cycle 
count ambiguity in the P OT b measurements, and there is also 
some uncertainty in the orbital period of IX Dra (see below) . 
However, the period uncertainty is in all cases negligible for 
the distance estimates given here. 



2.2.1 RX j '1715. 6+6856 

RX J1715. 6+6856 was not detected in 2MASS; we estimate 
Ks > 15.6 from the magnitude of the faintest stars detected 
in the field. At the period of RX J1715. 6+6856 (1.64 h), the 
donor sequence predicts Mk s = 8.74 for the secondary star 
(if this system is not a period bouncer; below we make the 
same assumption for SDSS J1730), giving d > 235 pc. 

Using the M v , diac -EW (H j3) relation, and EW(H(3) = 
62 A (from our spectroscopy) we obtain Mv = 10.2 ± 1.5 
(assuming that the disc dominates the flux in V), and with 
A v = 0.08 ± 0.04, we have d = 4001^0 P c - 



2.2.2 SDSS J 1730 

The discovery of SDSS J1730 was announced by 
ISzkodv et al. (2002), who observed a DN outburst, and mea- 
sured EW(Hf3) — 63 A (in quiescence) and an orbital period 
of roughly 1.95 h; the refined value is P orb = 1.837 + 0.002 h 
(John Thorstensen, private communication). As expected 
at this orbital period, SDSS J1730 is an SU UMa star 
(vsnet-alert 6798); the fractional superhump period excess 
indicates that it has not evolved to the period minimum 
dPatterson et al.ll2005l ). 

Kniggel ( 20061 ) gives Mk s =8.22 for a typical secondary 



at this period, and from the observed Ks = 15.22 ± 0.18 we 
therefore obtain a lower limit of d > 250 pc, and an estimate 
of d = 440l^u pc (using A v = 0.08 ± 0.04). 



3 The updated colour transformations are available at 



http://www.astro.caltech.edu/~jmc/2mass/v3/transformations/ 



4 The model presented in that paper that takes into account the 
spiral structure of the Galaxy. 
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The observed EW{Hj3) gives M v , dUc = 10.3 ± 1.5, and 
with V = 15.9±0.1, we have d — 130^60° pc (again assuming 
that the V-band flux originates only from the disc). Since 
this distance is only barely consistent with the robust lower 
limit, and since any white dwarf contribution to the V-band 
flux will imply that it is too low, we will adopt the donor- 
based distance estimate. 



2.2.3 RX J 183 1.1 +6511 

RX J1831. 7+6511 has P orb = 4.01 h (implying Mk s = 6.20 
for the secondary), K s = 14.96±0.14, and A v = 0.11±0.06 
(i.e. A Ks = 0.012 ± 0.007). This gives a lower limit of d > 
559 pc, and an estimate of d = 980lggQ pc. 

For comparison, we measure EW(H(3) — 14 A and 
V — 15.9 ± 0.2 from our 29 July 2006 spectra. Assuming 
the system was in quiescence, the Mv,disc~EW (H f3) rela- 
tion gives My < +6, and thus d > 900 pc. 



2.2.4 EX Dra 

Being a deeply eclipsing system, EX Dra has been the sub- 
ject of several det ailed studies. It is a U Gem star with 
Porb = 5.0385 h (|Fiedler et all Il997|). a secondary spec- 
tral type between M 1V and M2V (baptista et al.l l2000l ; 
iBillington et al.lll996l ; c.f. lHarrison et all 120051). and an or- 
bital inclination o f i = 85^ deg. ( Baptista et al.l I2OO0I ; 
iFiedler et al.l Il997l obtained a similar value). The out-of 
eclipse quiescent V magnitude is 15.3, and the system 
reaches V = 12.3 at maximum. At mid-eclipse in quiescence, 

V = 16.66 ± 0.07, I = 14.46 ± 0.05, and e ssentially all of 
the light originates from the secondary star dBaptista et al.l 
l200d ; iBaptista fc Catalanl l200ll ; IShafter fc Holland! 12003( 1 
Mid-eclipse photometry has provided distance estimates 
based on the photometric para l lax o f the secondary of 
240 ± 90 pc iShafter fc Holland! l2003h and 290 ± 80 pc 
(|Baptista et alj|2000h . " 

As a c heck, we can also use the CV donor sequence 
of iKniggel (l2006h to estimate the distance to EX Dra. 
We assume here that in quiescence, the secondary is the 
only source of mid-eclipse light. The sequence gives Mv = 
8.80 ± 0.37 for secondaries with spectral types in the range 
Ml to M2. This gives a distance of 360±g£ pc. At P orb = 
5.0385 h, predicted magnitudes for the secondary star are 
M v = 9.56 ± 0.96 and Mi = 7.51 ± 0.77 (the errors are 
derived from an error of ±1 spectral type in the sequence, 
and are probably conservative). If we use the measured pe- 
riod rather than spectral type, we find distance estimates 
of 260^90° an d 240lgQ° pc from the apparent mid-eclipse 

V and / magnitudes, respectively. The first of these esti- 
mates uses Av = 0.08 ± 0.04, while for the other two we set 
Av = 0.06 ± 0.03. The reason these distances differ is that 
EX Dra is n ot exactly on the donor sequence (see fig. 9 of 
iKniggd lipolj). 



with each other. We will adopt a distance of 26OI9Q P c f° r 
this system. 



2.2.5 IX Dra 



The M v 



-Porb relation gives d — 140I 50 pc. A likely 



reason for this small value is that the dependence on i given 
by equation [2] is too sensitive at large i, because radiation 
from the disc rim becomes significant (|Smakll 19941 ). 

Our three donor-bas ed di stance estimate s and those of 
IShafter fc Holland! (|2003l ) and IBaptista et al.l l|2000l ) are, of 
course, not independent, but they are all formally consistent 



iLiu et al. I <| 19991 ) confirmed IX Dra as a CV. It is a mem- 
ber of the group of SU UMa stars with unusually short 
outburst re currence times, som etimes referred to as ER 
UMa stars (|lshioka et alj l200ll ). The quiescent brightness 
is V = 17. 5, and maximum o f normal outburst is roughly 
V = 15.4. lOlech et all l|2004l ) measured a superhump pe- 
riod of 1.6072 h, and found a second photometric period of 
1.595 h in superoutburst light curves. If the second period is 
interpreted as P or t, the small fractional superhump period 
excess (e = (P sh — P or b) / Porb) implies a very lo w ratio of 
secondary to white dwarf mass (|Pattersonlll998j), and thus 
that IX Dra has a secondary of substellar mass joiech et all 
|2004| ). However, the period bouncer status of IX Dra is sus- 
pect for two reasons. The first is that the high outburst 
duty cycle (about 30%) implies a high M (however, note 
that DI UMa has very l ow e, and similar outburst charac- 
teristics; |Patterson| [l998l . More importantly, it relies on the 
1.595 h signal being the orbital period; a large amplitude (in 
this case about 0.05 mag) photometric orbital modulation is 
not usually seen during superoutburst, and the superhump 
modulation makes it dificult to measure other signals reli- 
ably. Confirmation of the orbital period is needed to decide 
whether this system is likely to be a period bouncer. 

The inclination of IX Dra is not known, but the fact that 
there is no eclipse (even at supermaximum) constrains it to 
< 70°. This gives a limit of d > 545 pc from the Mv,nax-P rb 
relation. 

2MASS provides only a lower limit on Ks and a 
marginal detection in H for IX Dra; we therefore use J = 
16.47 ± 0.12. If we used the period bouncer branch of the 
donor sequence, we would obtain a distance estimate that is 
inconsistent with the lower limit given above. We therefore 
use the pre-period minimum branch of the donor sequence 
(Mj = 9.85), despite the possibility that IX Dra is a period 
bouncer. In the direction of IX Dra, Av — 0.1 for distances 
greater than ~ 400 pc. Using the J-band offset of 1.56, we 
find d = 430i?9Q P c - 



3 THE SPACE DENSITY OF CVS 

Here, we give an estimate of the space density of CVs based 
on our sample (Section l3.1|) . and also consider, in Section l3T2l 
how large a population of systems with given (low) Lx could 
have escaped detection in the NEP survey. 

DNe typically have softer X-ray spectra in outburst 
than in quiescence, although a wide range in X -ray be- 
havio ur has been observed in different systems (e.g. IWarnerl 
1995). The RASS was conducted over a period of about 6 
months, and during this time the NEP area was observed 
many times. This means that, while the duration of the sur- 
vey was short compared to the outburst recurrence time- 
scale of faint DNe, high duty cycle, frequently outbursting 
CV, such as EX Dra, were probably observed both in out- 
burst and quiescence. We will not attempt to account for 
this complication. 
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Table 2. Orbital periods, X-ray fluxes, distances, and luminosities of the five CVs in the NEP survey area, together with the contribution 
each systems makes (if errors are ignored) to the mid-plane space density. Note that, although IX Dra is in the area covered by the 
survey, it is not included in the NEP sample. 



Object Porb/h Fjf/ergcm 2 s 1 d/pc L x /ergs 1 (l/V m ax)/pc 3 



EX Dra 


5.0385 


(6.4 + 0.7) x 10" 


14 


260±«° 


5 x 10 29 


8.8 x 10" 6 


RX J1831. 7+6511 
SDSS J1730 


4.01 
1.837 


(1.7 + 0.2) x 10" 
(3.4 + 0.3) x 10" 


13 
13 


980+^ 
U0±™ 


2 x 10 31 
8 x 10 s0 


1.4 x 10" 6 
3.1 x lO" 7 


RX J1715. 6+6856 


1.64 


(8.5 ± 1.7) x 10" 


14 


400+ 400 


2 x 10 30 


1.3 x 10" 6 


IX Dra 


1.595 


(2.2 ± 0.5) x 10" 


14 


190 


5 x 10 29 


5.4 x 10" 6 



We assume that the vertical density profile of CVs is 
exponential 



V) = D.- 



sin 6| ; 



[2 



[Xj + 2xj + 2) e 



(5) 



p(z) = p e~ 



t\/h 



(4) 



with z the perpendicular distance from the Galactic plane 
(i.e. z = dsin&, where b is Galactic latitude). The local space 
density is defined as po — p(0), the mid-plane value of p. 
We will ignore any radial dependence of p in Section 13.11 
(it certainly has a negligible effect on the results presented 
there), but include it in Section [3 . 2 1 (simply because it is not 
too computationally difficult). 

Widely used empirical values of the Ga l actic scale 
height of CVs, h, may be found in iPattersonl i| 19841 ) and 
iDuerbeckl l|l984T ). However, since all observational CV sam- 
ples contain a serious z-dependent bias, we will use theo- 
retical h values instead. Also, rather than assuming a single 
scale height for all CVs, we will follow IPretorius et al.l (|2007l ) 
in setting h = 120, 260, and 450 pc for long-period systems, 
short-period systems, and period bouncers, respectiveljjf]. 

CVs are expected to have X-ray emission that can be 
described as therma l bremsstrahlung with kT ro ughly be- 
tween 5 and 20 keV (jPatterson fc Ravmondll 19851 ); observa- 
tions indicate a larger variety of X-ray spectral shape s (e.g. 
IVrtilek et af1ll994l ; lRlchmanlll996l ; [Navlor et al.lll988l ). The 
range in Fx resulting from varying kT from 5 to 20 keV for 
a thermal bremsstrahlung spectrum is about as large as the 
observational error (which is an insignificant contribution 
to our uncertainty in po). We will therefore simply assume a 
kT — 10 keV thermal bremss trahlung spectrum for all CVs. 
We follow iHenrv et all (120061 ) in quoting X-ray fluxes and lu- 
minosities in the 0.5-2.0 keV band, calculated from ROSAT 
PSPC count rates in the 0.1-2.4 keV band. 



3.1 The l/Vmax method 

The space density is found by counting the CVs inside the 
volume observed, while accounting for both the fact that our 
sample is flux- rather than volume-limited, and the variation 
of space density with z. The method involves essentially al- 
lowing the 'volume limit' to vary according to the luminosi- 
ties of the systems in our sample, while the dependence on 
z is taken care of by defining a 'generalised' volume 



5 This is clearly still a very simplified picture, but it is an im- 
provement over assuming a single small scale height for the entire 
CV population. 



(e.g. ISchmidtJ Il968l : iFeltenl ll976J : IStobie et al l Il989l : 
iTinnev et all 1 19931 ) . where Q. is the solid angle covered by 
the survey and Xj = dj \ sinb\/h, with dj the maximum dis- 
tance at which a CV included in the sample could have been 
detected (given its Lx and the survey flux limit). Equation[5] 
assumes that p is a function only of z, and that this depen- 
dence is given by equation^ The index j represents the CVs 
in our sample. Because both b and the flux limit are variable 
over fi, we compute each Vj as a sum over solid angles 8£l 
small enough so that b and the flux limit (and thus dj) can 
be assumed to be constant over each SQ. This is done by 
dividing the NEP area into a 36 x 36 pixel grid, po is then 
the sum of the space densities represented by each CV, i.e. 



po 



The l/Vmax values for the 5 CVs, obtained by using the best 
distance and Fx estimates, are given in Table [2] together 
with other relevant parameters. 

We neglect extinction in finding the maximum distance 
at which a source could be detected from its estimated dis- 
tance (for which extinction was included) and observed Fx ■ 
This is reasonable because extinction is low in the NEP area 
(the absorbed Fx is at least 0.78 times the unabsorbed flux 
for all positions in the survey volume), and because, for the 
distances of most of our sources, the column densities are 
already at close to the total Galactic values. In other words, 
at distances beyond our sources, extinction is approximately 
constant, so that the maximum volumes we obtain are un- 
affected if extinction is ignored. 

Both distance errors and the small sample size con- 
tribute significantly to the overall uncertainty on po, al- 
though the latter effect dominates. In order to account fully 
for all sources of statistical error (including the uncertainty 
in distances, and also the error on Fx), we determine the 
probability distribution of po via a Monte Carlo simulation. 
Thus we create a large number of mock CV samples with 
properties designed to fairly sample the full parameter space 
permitted by the data. We then compute an estimate of po 
for each mock sample, thus gradually generating a probabil- 
ity distribution for po. 

The key step in this Monte Carlo simulation is the gen- 
eration of appropriate mock samples. This is done by treat- 
ing each CV in the real sample as a representative of a larger 
population of similar CVs, and allowing each of these pop- 
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ulations to contribute to each mock sample. In creating a 
particular mock sample, we therefore consider each CV in 
the real sample in turn, and generate a mock counterpart to 
it by drawing the X-ray flux and distance from the appro- 
priate probability distributions. This takes care of the ob- 
servational uncertainties in Fx, as well as the uncertainties 
in distances. However, there is also Poisson uncertainty as- 
sociated with the small sample size, which would be ignored 
if each mock sample contained exactly the same number of 
CVs as the real sample. We deal with this by assigning a 
weighting factor to each counterpart CV in each mock sam- 
ple. The weighting factor needs to be drawn independently 
for each counterpart CV from the probability distribution of 
the expected number of sources belonging to this population 
in a NEP-like survey. This probability distribution can be 
constructed from Bayes' theorem. If the expected number 
of sources in a NEP-like survey is p; for a particular popula- 
tiorfl the actual number observed in a particular survey will 
be N i,s, where this is drawn from the Poisson distribution 
P(N i, s \[i). Then Bayes' theorem gives 

P(p\N obs ) oc P(AU»P( M ), 

where P(p) is the prior on p. Note that, in practice, N a t s = 1 
by construction for each population represented in the real 
sample (since each CV is taken to represent a population 
of similar CVs). Also, since /i is a scale factor for the space 
density of each population, the appropriate uninformative 
prior to use is P(p) = l/fi (e.g. |Jeffrevsl[l9 6lf1. 

To summarise, we draw a random distance and Fx from 
the appropriate probability distribution functions (assuming 
that the error in Fx is Gaussian) for each of the systems 
in the sample, compute Vj, generate fj,j, and obtain po by 
summing pj /Vj over j. Repeating this a large number of 
times produces a probability distribution for po. Fig.|5]shows 
the po histogram obtained from this calculation. The mode, 
median, and mean of this distribution are 4.7 x 10 -6 , 1.1 x 
10 -5 , and 2.4 x 10 _5 pc -3 , and are marked by solid lines. 
The dotted lines show the symmetric 1-a confidence interval 
around the median, which is 4.4 x 10~ 6 to 3.5 x 10~ 5 pc -3 . 

A reason for some concern is that IX Dra was only just 
too faint to be included in the NEP sample (the fact that 
it might be a period bouncer adds to this problem). More 
importantly, a single system, EX Dra, dominates the space 
density estimate. Our choice of scale heights is also non- 
standard. To investigate the sensitivity of the result to these 

6 Note that p does not need to be an integer. 

7 We carried out several tests to check that this technique pro- 
duces reliable error estimates. In the tests, we take a space den- 
sity, luminosity function, and flux limit as input, and use these to 
simulate observed samples. The po distribution implied by each 
sample is then calculated in the same way as above. We found that 
for ~68% of these po distributions the true input po is contained 
in any particular 1-a interval. This holds when the expected num- 
ber of detected sources is as small as 4, but also for much larger 
samples, when the distributions become narrower. We ignored 
Galactic structure in these tests (i.e. p was taken to be constant), 
but experimented with several different (simple) luminosity func- 
tions, as well as with including errors in luminosity. We also ran 
some simulations to verify that a uniform prior is not appropri- 
ate in this calculation. These numerical tests indicate that the 
method we use is reliable, and that we have adopted a suitable 
prior. 



3 




lO" 6 SxlO" 5 3xl0" 5 4xl0" 5 

Po/P c " 3 



Figure 6. The po distribution resulting from our simula- 
tion. This sample includes RX J1715. 6+6856, SDSS J1730, RX 
J1831. 7+6511, and EX Dra; we used scale heights of 260 pc and 
120 pc for short- and long-period systems, respectively. Solid lines 
mark the mode, median, and mean at 4.7 x 10" 6 , 1.1 x 10" 5 , 
and 2.4 X 10 — 5 pc - 3 , and dotted lines show a 1-a interval from 
4.4 x 10" 6 to 3.5 x 10" 5 pc~ 3 . 



Table 3. Systems are identified as 1 — EX Dra, 2 — RX 
J1831. 7+6511, 3— SDSS J1730, 4— RX J1715. 6+6856, and 5— IX 
Dra in the first column, hi and h 2 are the scale heights, in parsecs, 
assumed for short- and long-period systems, respectively. The val- 
ues given for po are the medians of the po distributions, and the 
errors represent symmetric 1-a confidence intervals, around the 
median. 
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systems 1, 2, 3, 


4; 


hx = 260, h 2 = 120 


1.1" 


^2.3 
-0.7 


X 10" 


5 


systems 1, 2, 3, 


4; 


hi = h 2 = 120 


1.5; 


H2.5 
-0.8 


x 10" 


5 


systems 1, 2, 3, 


4, 


5; /ii = 260, h 2 = 120 


2-i; 


1-4.2 
-1.3 


x 10" 


5 


systems 2, 3, 4; 


hi 


= 260, h 2 = 120 


3.1; 


H5.5 
-1.9 


x 10" 


6 


systems 1, 2; h 2 




120 (long-period) 


7.7; 


US 
-5.1 


x 10" 


6 


systems 3, 4; h 2 




260 (short-period) 


i.i" 


15.2 
-1.0 


x 10" 


6 



factors, we list in Table |3] the space densities that would be 
obtained if we (i) used a single scale height of 120 pc, (ii) 
included IX Dra, and (iii) excluded EX Dra. Note that as 
different as the best estimates (we give the median of each 
distribution) are, they are all consistent with each other. 
Table [3] also gives the space densities for short- and long- 
period systems separately. 



3.2 Upper limits on the space density of an 
undetected population 

An unavoidable assumption of the 1/Vmax method is that 
the detected objects are representative of the luminosity 
function of the true underlying CV population. No evidence 
of even a large population of CVs is expected to show up 
in a flux-limited survey, if such a population is sufficiently 
faint. 

The faintest meas ured values of Lx for CVs are a few 
times 10 29 ergs _1 s fe.g. lVerbunt et al ]|l997l ; IWheatlev" et alj 
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2000; ISchwope" et al . 20070 

However, we do not really know 
how faint CVs can be: the faintest expected secular lumi- 
nosities can be estimated from the gravitational radiation 
M, but, since faint CVs are DNe, they spend most of their 
tim e at fainter l umino sities . 

iRiaz et all i|2006j ) find lO^ergs" 1 < L x < 10 33 ergs -1 
for a large sample of nearby dwarfs with spectral types be- 
tween K5 and M6; brown dwarfs have Lx < 10 29 ergs _1 
(often several orde rs of magnitude less than 10 29 erg s _1 ; e.g. 
IStelzer et al]l2006h . CV secondaries have shorter spin peri- 
ods than expected for the objects in these studies, but, even 
so, we cannot be confident that the secondary star alone 
guarantees Lx > 10 29 ergs -1 . This is especially true in the 
case of period bouncers. 

We use another Monte Carlo simulation to find the 
upper limit on the space density of a hidden population, 
assuming that all members of this population have the 
same (low) Lx- The Galaxy is modelled as an axisymmet- 
ric disc, with exponential vertical and radial (r) density 
profiles. The galactic c entre distance is taken as 7 620pc 
l|Eisenhauer et all [2005) . and the small offset of the Sun 
from the galactic plane is neglected. At least part of a hid- 
den population may consist of pre-period bounce systems, 
but we conservatively assume a single scale height of 450 pc 
(appropriate to period bouncers). The radial scale length is 
taken as 3 000 pc. The local space density is then po, where 
p = po exp(— 1«|/450 pc), for r near 7620 pc. 

Because our upper limits will be systematically low if 
ex tinction is neglected, we include it here, using the model 
of [A mores fc Lepinel (|2005h to find Nh for every simulated 
system. 

We then find the value of po for which the predicted 
number of faint systems detected in the NEP survey is 3 
(detecting such systems is then a 2-<r result). Fig. [7] shows 
the maximum allowed po as a function of L x for CVs that 
make up the hidden population. Results from the simulation 
are plotted as a fine histogram, while the bold curve is a fit 
to the data. The fit is given by 



p max = 5.7 x lo-^wnrVgs- 1 )- 



pc 



Thus for L x = 2 x 10 29 ergs -1 , we have p < 2 x 10 -4 pc" 3 . 
Note that here po refers only to the hypothesized undetected 
population, and does not include the contribution from the 
brighter, detected systems. 



4 DISCUSSION 



The po measurement from the NEP survey is consistent with 
several previous observational estima tes, e.g. ~ 1 0~ 5 pc~ 3 
l|Pattersor]|l99sl ) and ~ 3 x 10" 6 pc" 3 (IWarneifeOOlh A par - 
ticularly relevant study to compare to is iHertz et all l|l990l ). 
since it was also based on a sample of 4 X-ray selected CVs. 



8 IVerbunt et ail lll997f) list L x = 3.16 X lO^ergs" 1 for AY 
Lyr, but this resulte d from a very low distance estimate of 57 pc 
llSproats et al.|[l996h . The distance is based on a measurement of 
K = 13.08, taken during outburst. 2MA SS phot o metry indicates 
Kg > 15.6 in quiescence, and using the iKniggel d20odh sequence 
(with no offset), we find d > 277 pc, implying that the system is 
not particularly faint in X-rays. 




39.8 29.4 
lo g(Lx/ er g s ~') 

Figure 7. The upper limit on the mid-plane space density as a 
function of X-ray luminosity for an unseen population of CVs. 
The data from the simulation are shown as a fine histogram, and 
a fit is over-plotted as a bold curve. 



Hertz et al. find po ~ (2 — 3) x 10 _5 pc -3 (where the un- 
certainty includes only the possible incompleteness of the 
sample), again consistent with the estimate presented here. 

The most important concern with our measurement is 
that the NEP sample is so small that it may not fairly rep- 
resent the underlying CV population. This is apparent from 
the fact that we estimate a higher space density for long- 
than for short-period systems. The reason is that our esti- 
mate is dominated by the long-period system EX Dra, which 
has an unusually low Lx for its orbital period, probably be- 
cause of its very high inclination. Note, however, that if the 
X-ray luminosity of EX Dra had been higher, that would 
only reduce po. 

Theoretical predictions of the CV space density are typ- 
ically higher than most observational values. A space den- 
sity of (0.5 - 2) x 1 0~ 4 pc~ 3 results from the CV birth 
rate of Ide Kool (|l992T ). together with the assumptions that 
the age of the Galaxy is f0 Gy, and that the lifetime 
of a CV is longer than t hat J Warner! [1991 ). When the 
model population of iKolbl (|l993 ) is normalised to the for- 
mation rate of single white dwarfs, the predicted mid- 
plane CV space density is 1.8 x 10 _4 pc -3 . We cannot 
rule out these values, but our data do show that if po is 

pc 

ergs 

pected that the mass transfer rates of CVs fluctuate around 
the secular mean on time-scales muc h longer than the to- 
tal observational base line (see e.g. [H ameury ct al. 1989; 
Wu et al.lll995l;lKing et al.ll 19951; [McCormick fc Franklll998l ; 
Biining fc Ritterj|2004l ; lRitter et al.ll2000l ') . Because the time- 
scale of these M cycles is not accessible to observations, we 
can only note that surveys useful for finding CVs typically 
have no sensitivity to objects that are not currently accret- 
ing. 

The fact that stellar populations of different ages do 
not have the same Galactic distribution has to be kept in 
mind when interpreting (and, indeed, measuring) space den- 
sities. While the situation is clearly not really as simple 
as three distinct sub-populations with three scale heights, 
long-period CVs make up a significantly younger popu- 



really as high as 2 x 10 4 

population must have Lx < 2 x 10 29 ergs _1 . It is sus- 
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lation than, say, period bouncers. One should therefore 
clearly not expect the ratio of the mid-plane space den- 
sities of these sub-populations to be the same as the ra- 
tio of their total Galactic numbers. This may be a triv- 
ial point, but it has been completely ignored in the liter- 
ature. If, e.g., we assume that 70% of all CVs in the Galaxy 
are period bouncers, and only 1% are long-period systems, 
and if the space densities of these two sub-populations are 
given by p pb {z) = p pb {0) exp(-|z|/450 pc) and pi p (z) = 
pi p (0) exp(— 2; |/120 pc), respectively, then (from integrating 
the vertical density profiles and assuming the radial and az- 
imuthal dependences are the same in both cases) we have 
p pb (0) = 70(120/450)p ip (0) ~ 20p ip (0). We can also use this 
argument to place a tentative limit on the ratio of the total 
numbers of long-period CVs and periods bouncers. If we as- 
sume p pb (0) < 2 x 10~ 4 pc~ 3 and p tp (0) > 2.0 x 10~ 6 pc~ 3 , 
then this ratio must be less than about 30. 

There are several ongoing surveys that will produce 
large new samples of CVs, e.g. the Hamburg Quasar Survey, 
the Sloan Digital Sky Survey, and the INT Photometric Ha 

Gansicke et al.l 



Surv e y of the Northern Galactic Plane (e . L. t 
20021 ; ISzkodv et ail |2002| ; IWitham etaLl 120061 ). Although 



these samples will all be subject to more complicated selec- 
tion effects than the NEP sample, they will be much larger, 
which will remove the main shortcoming of the present study 
(although it should be kept in mind that space density esti- 
mates based on l arger samples can st ill be dominated by 1 
or 2 systems; e.g. ISchwope et al. 2002). The Palomar Green 
Survey CV sample (e.g. |Ringwald 1993) is already available, 
and is also suitable for a rigorous measurement of po, since 
it is complete and has well-defined selection criteria. 



5 CONCLUSIONS 

Using the complete X-ray flux-limited CV sample detected 
in the ROSAT NEP survey, we have estimated the space 
density of CVs bright enough to be represented in the survey, 
and placed limits on the size of a fainter population that may 
have remained undetected. 

We have presented observations of two newly discovered 
CVs, namely RX J1715.6+6856 and RX J1831. 7+6511, and 
measured orbital periods of 1.64 h and 4.01 h respectively 
for these systems. RX J1715. 6+6856 and RX J1831. 7+6511 
are probably both DNe. This completes the identification of 
CVs in the ROSAT NEP survey. 

Distances were estimated for the 4 systems in the NEP 
sample, as well as for IX Dra, which is in the area covered 
by the survey, but below the survey flux limit. 

With the assumption that the 4 CVs included in the 
NEP survey are representative of the intrinsic CV popu- 
lation, the space density of CVs is l.l^n'? x 10 -5 pc -3 . 



Space densities as high as the theoretically predicted 



2 x 10 4 pc 3 are excluded by our data, unless the dom- 



inant CV population is fainter than ~ 2 x 10 ergs - in 
X-rays. 

We also emphasise the fact that long-period CVs, short- 
period CVs, and period bouncers have different Galactic 
scale heights; this needs to be accounted for when mid-plane 
space densities are translated to total Galactic number of 
systems. 
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